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ABSTRACT 

jy\ ' The assumption of Gaussianity of the primordial perturbations plays an im- 

portant role in modern cosmology. The most direct test of this hypothesis 
consists in testing Gaussianity of the CMB maps. Counting the pixels with 
the temperatures in given ranges and thus estimating the one point proba- 

m ■ 

t^^ ■ bility function of the field is the simplest of all the tests. Other usually more 

o ■ 

complex tests of Gaussianity generally use a great deal of the information 
already contained in the probability function. However, the most interesting 
O , . outcome of such a test would be the signal of non-Gaussianity independent 

Q ' of the probability function. It is shown that the independent information has 

H ■ 

"^ ' purely morphological character i.e. it depends on the geometry and topology 

of the level contours only. As an example we discuss in detail the quadratic 
model V = u + Ci{\i} — 1) (u is a Gaussian field with u = and < iP' >= 1, 

\^ ' a is a parameter) which may arise in slow-roll or two-field inflation models. 



We show that in the limit of small amplitude a the full information about 
the non-Gaussianity is contained in the probability function. If other tests are 
performed on this model they simply recycle the same information. A simple 
procedure allowing to assess the sensitivity of any statistics to the morphologi- 
cal information is suggested. We provide an analytic estimate of the statistical 
limit for detecting the quadratic non-Gaussianity olc as a function of the map 
size in the ideal situation when the scale of the field is resolved. This esti- 
mate is in a good agreement with the results of the Monte Carlo simulations 
of 256^ and 1024'^ maps. The effect of resolution on the detection quadratic 
non-Gaussianity is also briefly discussed. 

Key words: Cosmic Microwave Background Radiation - Cosmology. 
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2 S.Shandarin 

1 INTRODUCTION 

The quest for the physical mechanism of the generation of the initial inhomogeneities along 
with the measurements of the major cosmological parameters {Hq, Qa, Qcdm, ^b, P{k), etc) 
is one of the most important problems in modern cosmology. The standard inflationary 
model predicts the primordial fluctuations were Gaussian random fields ( |Guth fc Pi 1982| ; 



Hawking 1982| ; |Starobinski 1982] ; [Bardeen, Steinhardt &: Turner 1983| ). In agreement with 



the theory the current observations provide little evidence for deviations from Gaussianity. 
The majority of the tests of Gaussianity in the COBE maps ( |Golley, Gott fc Park 1996| ; 
Kogut et al. 1996|; [berreira, Magueijo fc Gorski 1998|; [Novikov, l:'eldman fc Shandarin 1999|; 



Pando, Valls-Gabaud fc Fang 1998| ; [Bromley fc Tegmarkl999| ; [Banday, Zaroubi fc Gorski 



20001 ; [Mukherjee, Hobson fc Lasenby 2000| ; [Barreiro et al. 2000| ; |Aghanim, Forni fc Bouchet 



2001 ; Phillips fc Kogut 2001| ) have resulted in the general agreement that all non-Gaussian 



signals were of noncosmological origin Q. This was not perhaps a surprise because of a very 
large physical scale corresponding to the COBE resolution (^ 7 deg). Recent studies of the 
AT/T maps on a degree and sub-degree scales also showed no significant deviations from 
Gaussianity ( [Park et al. 200 1| ; fShandarin et al. 200 1| ; |Wu et al. 20011) . 



Nevertheless, the question of possible non-Gaussianity in the CMB maps is very impor- 
tant because of the following reasons. First, a detection of a non-Gaussian component in the 
primordial fluctuations may profoundly affect modern cosmology ruling out some models 
of the early universe and boosting the others (see e.g. Turner ( |1997| ); Vilenkin & Shellard 
( p.994| )). Second, Gaussianity is a key underlying assumption of all experimental power spec- 
trum analyses to date, entering into the computation of error bars ( [I'egmark 1997| ; Bond fc 



J arte 199^ ), and therefore needs to be observationally tested. In addition, the hypothesis of 



the Gaussianity of the initial perturbations enters in the studies of the large-scale structure, 
clusters of galaxies, formation of galaxies, Lya forest, etc. Finally, studying Gaussianity of 
CMB maps may reveal otherwise undetected foreground contamination. 

Many tests for non-Gaussianity have been suggested. For instance, in the recent paper 
by Wu et al. ( p001|) the authors employed a total 82 hypothesis tests for Gaussianity, 
although, as the authors noted, the tests were not independent. Obviously the question 
arises how independent the different tests were. A related issue concerns the possibility of 
constructing a set of independent tests. The full solution to this problem is beyond the scope 



* However, Magueijo (20001) still has a 97% confidence level that the signal is not due to systeniatics. 
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Morphological non-Gaussianity 3 

of this paper but we describe a simple procedure that allows to assess the independency of a 
test from the probability function. The probability function, f{u) also called the probability 
density function tells what is the probability that a randomly chosen point of the field, u{x, y) 
has a certain value. In practice f{u)du is often estimated by counting the cells having the 
values of the field between u and u + du. Gaussian fields have the Gaussian probability 
function. The choice of the probability function (PF) as the reference statistic is due to both 
its conceptual and practical simplicity. Generally speaking the more sensitive a test to the 
statistical information different from the probability function the more useful it is. 

The hierarchy of n-point correlation functions or equivalently n-spectra is also widely 
used for testing non-Gaussianity. However, compared to PF this is not as easy task from 
both conceptual and practical points of view. One reason is the multidimensional character 
of the n-point functions. The n-point function as well as n-spectrum in the two-dimensional 
space is a function of 2r2 — 3 variables. Thus, testing the lowest order functions sensitive 
to non-Gaussianity, three-point correlation function or bispectrum, one has to deal with a 
function of three variables. So far pragmatic solutions to this problem were either computing 
a small number of particular cuts in the three-dimensional space or introducing some integral 
quantities. Both shortcuts obviously result in incompleteness of the test. The other reason 
is purely computational: computing of a n-point function on the grid with Npi^ pixels using 
current methods requires 0{Np) operations which is already prohibitive for current fairly 
small maps (COBE, QMASK, Maxima I) even for ra = 4 — 5. Although, a clever technique 
using kd-trees can potentially reduce it to 0{Npix{logNpix)"'~^) ( [Szapudi et al. 20001) it has 
to be developed yet. 

The n-point functions carry information about the maps in highly redundant and diluted 
form. In order to see this let us consider a large two-dimensional map obtained observation- 
ally or from a theoretical simulation. Obviously all information about the map can be stored 
in the form of a function of two variables (e.g. the map itself). The two-point function of the 
map is a function of one variable and thus considerably reduces the information about the 
map by loosing the phase information. In general case the three-point function also consid- 
erably reduces the information about the map but in contrast to the two-point function it 
increases the dimensions of the space that means a huge dilution of information. The four- 
point function in the two dimensional space is a function of five variables meaning that it 
dilutes the information even more than the three-point function. In general case, the higher 
order of the n-point function the more diluted is the information about the field. 

© 0000 RAS, MNRAS 000, 000-000 



4 S.Shandarin 

It is well known that every n-point function affects the probability function, see e.g. 
White (1979) and Balian and Schaeffer (1989). It means that the PF carries some information 
about every n-point function. Reversing this statement one can say that no n-point function 
carries information that is completely independent from the PF. Thus, one may also ask 
what information is stored in n-point functions which is independent of the PF and whether 
it is possible to extract it or at least to assess it. Obviously, the same question must be 
addressed not to only the n-point functions but to all other statistics. These issues are 
discussed bellow. 

The PF or equivalently the cumulative probability function (CPF) [y is not only the 
simplest conceptually but also most efficient numerically. Computing this statistic requires 
only 0{Npix) operations. The only problem is that the Gaussian PF does not guarantee the 
Gaussianity of the field. Therefore, some additional statistical information is badly needed 
in case when the PF of the field is Gaussian since if the PF is non-Gaussian the non- 
Gaussianity is already detected. The next step obviously would be the identification of the 
physical process responsible for the non-Gaussianity but ffist it must be detected. Thus, if 
the PF is Gaussian the additional information must be independent of that contained in the 
PF. We will show that such information has purely morphological character. This means that 
it is completely determined by the geometric and topological statistic of the excursion sets. 
Thus, a set of morphological parameters based on Minkowski functionals becomes a natural 
choice of the statistics that is sensitive to non-Gaussianity and completely independent of 
the PF provided that proper parameterization is used. 

A particular kind of non-Gaussianity known as the quadratic model has recently attracted 
much attention ( [Coles fc Barrow 1987| ; [Luo fc Schramm 1993| ; [Matarrese, Verde fc Jimenez 



2000t [Verde et al. 200q ; P/erde 200 1[ ; [Verde et al. 2001|) . One reason is that it could 



be generated by plausible physical mechanisms in the early universe ([Falk, Rangarajan fc 



Srednicki 1993[ ; [Gangui et al. 1994[; [Luo 199^ ). The other is a relative ease of its analysis. In 



this paper we show that the simplest test for Gaussianity, the probability function, provides 
also the complete statistical information in the most interesting case of small amplitudes. It 
means that other tests if applied to this model at best only recycle a part (probably small) 

T Here we would like to emphasize the informational content of a non-Gaussian statistic and therefore do not distinguish the 
PF and CPF assuming that both contain the same information. 
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Morphological non-Gaussianity 5 

of this non-Gaussian information. In the general case of arbitrary amphtude the set of global 
Minkowski functionals completely characterize the statistical properties of this field. 

The rest of the paper is organized as follows. We describe the set of morphological quan- 
tities in Sec. 2. Sec. 3 describes a particular parameterization of the morphological statistics 
that makes them PF- independent. As an illustration we discuss a couple of simple non- 
Gaussian models one of which is the quadratic model that is often used in cosmology in 
Sec. 4. Then, in Sec. 5 we describe a class of the simplest non-Gaussian fields which can 
be called trivial non-Gaussian fields. Section 6 describes a few simple estimators of the 
amplitude of the quadratic non-Gaussianity. We describe simple Monte Carlo simulations 
modeling the detection of the quadratic non-Gaussianity in Sec. 7. Finally, we discuss the 
results in Sec. 8. 

2 MORPHOLOGICAL QUANTITIES IN TWO DIMENSIONS 

Morphology of two-dimensional random fields can be conveniently described in terms of 
geometric and topological properties of the regions bounded by the contours of constant 
level. There is a particularly useful set of quantities called Minkowski functionals ( [Minkowski 



which have very simple geometric and topological interpretations. For each isolated 



region bounded by a contour there are only three scalar Minkowski functionals: the area 
within the boundary, a^, its perimeter or the contour length, q, and the Euler characteristic 
or genus, gi which is 1 — nhi where rthi is the number of holes in the region. 

Minkowski functionals are additive quantities and therefore can be easily calculated for 
any set of regions if they are known for each region. In particular, the global Minkowski 
functionals, i.e. the total area. A, contour length, C and genus, G of the excursion set: 

A = ^a, , C = ^c, , G = Y.9^ (1) 

i i i 

are often used (|Gott et al. 1990 ; Winitzki fc Kosowsky 1997 ; Schmalzing fc Gorski 1998 ; 



Novikov, Feldman fc Shandarin 1999| ; [Novikov, Schmalzing fc Mukhanov 2000] ; [Wu et 



al. 20011 ; [Shandarin et al. 200 1|) . 



The total area of the excursion set, A is obviously the CPF of the field: A{u) = F{u) = 
IT f{u')du'. The Euler characteristic or genus have been used in cosmology for a number of 
years ( Poroshkevich 1970| ; |Gott et al. 1986t [Coles 1988|) . 



The first time the set of global Minkowski functionals was introduced into cosmology with 
the reference to their significance in differential and integral geometry by Mecke, Buchert 
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6 S.Shandarin 

& Wagner (| 1994|) and Schmalzing & Buchert ( | 1997|) . In particular, they emphasized a 



powerful theorem by Hadwiger (| 1957|) stating that under rather broad restrictions the set 
of Minkowski functionals provides a complete description of the morphology (for further 
discussion see e.g. Kerscher (| 1999|) ). 

In addition, the Minkowski functionals of the largest (by area) region {Ap, Cp, and Gp) 
give accurate description of the percolation phase transition ( |Yess fc Shandarin 1996| ). At 
percolation the regions merge into one region that spans throughout the whole space of the 
field. Percolation phase transition is sensitive to some types on non-Gaussianity ( [Zerdovi"cE 



r982t [Shandarin 1983| ; [Shandarin fc Zel'dovich T983| ; [Shandarin fc ZerdovichT98^ ). 



3 PARAMETERIZATION 

Often the level of the excursion set, u is used to parameterize Minkowski functionals. The 
global Minkowski functionals of a Gaussian field as functions of the level are given by simple 
analytic formulae ([Longuet-Higgins 1957[ ; pott et al. 19901) 



A(u) = 5[l-erf(^^ 

where R = ^/2/ai is the scale of the field; ai is the rms of the first derivatives (in statistically 
isotropic fields both derivatives du/dx and du/dy have equal rms). It is assumed that the 
field is normalized: < u >= and < v^ >= 1. The Minkowski functionals as functions of 
the level are shown in the left hand side panels of Fig.l. 

The parameterization by the level is useful for many applications. However, for the study 
of the morphology of the fields and Gaussianity in particular the parameterization by A is 
much better because it makes the statistics independent of the PF and considerably less 
correlated with each other ( phandarin et al. 2001[) . 

The total area of the excursion set A of the field was suggested to parameterize other 
quantities ( [Yess fc Shandarin 1996| ; [Shandarin et al. 2001[ ). Parameter A is a single valued 



function of z/, the parameter used in the most papers studied genus (Park et al. ( [ 2001 



and references therein) and therefore every function of u can be also expressed as a function 
of A. However A is more directly related to the morphology of random fields. In addition, 
being equal to the cumulative probability function (CPF) it has a very simple statistical 
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Morphological non- Gaussianity 7 

and geometrical meaning. Imagine that the excursion set {u > Uc) is painted black while the 
rest of the map {u < Uc) remains white. Then the fraction of the area in black equals A. 
The right hand side panels in Fig.l show the level, u, total contour length, C and genus, G 
as a function of the total area of the excursion set A for a Gaussian field. 

It has been noticed that percolation properties of the field can be useful for detecting 
some types of non- Gaussianity ( |Zerdovich 1982| ; ^handarin 1983| ; |Shandarin fc Zel'dovich 



1983| ; |Shandarin fc Zel'dovich 1984| ; [Klypin fc Shandarin 1993]) . The Minkowski functionals of 



the largest by area region are excellent parameters to characterize the percolation properties 
of the field. In contrast to the global Minkowski functionals they are not known in an analytic 
form. Figure 2 shows the Minkowski functionals of the largest region for a Gaussian field 
obtained from a large ensemble {N = 1000) of large Monte Carlo simulations (1024^ grid). 
It should be stressed that the major reason of using A as an independent parameter 
consists in isolating independent morphological information that is not present in the prob- 
ability function (PF). For example, C = C{A) carries only PF-independent information, 
while C = C{u) mixes it up with the information stored in the PF. 

4 EXAMPLES OF NON-GAUSSIAN FIELDS 

It is interesting to compare how some of the n-point functions and morphological charac- 
teristics signal the presence of non-Gaussianity in the field. We consider two examples with 
quadratic and cubic non-Gaussianity. The former has been suggested as a model having 
plausible physical mechanisms producing small deviations from Gaussianity ([balk, Kangara- 



Ijan fc Srednicki 1993| ; [Gangui et al. 19U^ ; [Luo 1994] ) . The latter has no physical motivations 



and is taken as a toy model only. 

4.1 Cubic Model 

First we consider a transformation 

v = u + (3u^, (3) 

where f3 is assumed to be positive, which guarantees monotonicity of the mapping. The 
parent field u is assumed to be a normalized Gaussian field with < n >= and < u^ >= 1. 
The CPF of the resulting field is 

Fiv) = Fcius), (4) 
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8 S.Shandarin 

where Ug is the solution of the cubic equation Ug + /3ul = f q Ahhough the solution could be 
obtained in a closed form the expression is quite cumbersome therefore we take the limiting 
case of small /? 

Thus, the CPF becomes 

/^„3 \ 1 ^ 

(6) 



F{v) = - 1 1 - erf 



1 / (3v^ 

V 



v/2 V l + 3(3v\ 
and the PF can be obtained by differentiation: f{v) = —dF/dv. 

The two- and three-point correlation functions of this field can be easily obtained for an 
arbitrary (3 

6(1,2) = (1 + 6/3)^,(1, 2) + 3/3^^,(1, 2) [3 + 2^2(^^2) 
^3(1, 2, 3) = (7) 

where ^(1, 2) is a shortcut to ^(ri, r2). The function ^g is the two-point correlation function 
of the parent field u. To the linear order in (3 the two point function remain the same by 
the form but acquires a different normalization (1 + Q(3). The reason why the three point 
function does not show the presence of non-Gaussianity is the symmetry of the mapping: 
n — ^ f is mapped by an odd function. 

Global Minkowski functionals C{A) and G{A) as well as the percolation statistics Ap = 
Ap{A), Cp = Cp{A) and Gp = Gp{A) remain exactly same as in the parent Gaussian field 
since the transformation simply relabels the levels without changing the contour lines. How- 
ever, if they were expressed as functions of the level {G{v) and G{v)) they would differ 
from Cg and Gg because some non-Gaussian signal leaks into them from the CPF. We will 
illustrate this point when discuss the quadratic model. 

The full information about non-Gaussianity of the field is obviously stored in one number 
(the value of f3) and can be obtained from one point statistics (the PF or equivalently CPF) 
however the non-Gaussianity is not detected by the odd order n-point statistics. In a generic 
case of monotonic mapping v = 0(m) where is a monotonic but not necessarily odd function 
of u all n-point functions would detect non-Gaussianity. 

■f For the positive /3 there is only one real solution Ug. 
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Morphological non- Gaussianity 9 
4.2 Quadratic Model 

As it was mentioned before the quadratic non- Gaussianity is particularly popular in cosmol- 
ogy. The quadratic model is the sum of a Gaussian field and its square 

V = u + a{u'^ -I). (8) 

In this paper we assume that u is normalized to unity: < u >= and < u"^ >= 1. If one likes 
a different normalization < u' >= but < u'"^ >= a'"^ and the quadratic transformation in 
the form 

v' = u' + a'iu'^ - a'2) (9) 

then the relations between the parameters of the two transformations are as follows 



u 



a M, 



a f , 



a 



a/a'. 



For certainty, without losing generality we will assume a > 0. Solving eq. 
denoting the solutions as ui and U2 {u2 < mi) one obtains 

'l + Aa{a + v) - 1 



Ml 



1 +4a(a + f) - 1 



U2 



2^a ' 2V2a 

The CPF of the field can be written then as 



1 



'I +Aa{a + v) - 1 



l + Aa{a + v) - 1 



(10) 
for u and 

(11) 

(12) 



F{v) = 1 - - erf ^ -^ - erf , 

^ ^ '2[ \ 2V2a J \ 2V2a 

Differentiating it with respect to v one easily obtains the PF (P^uo fc Schramm 1993 ; Matar 

rese, Verde fc Jimenez 200"0| ) ^ 



f{v 



1 exp 



8«2 ) -r exp 



(-^l+4a(a+i>)-l)^ 
8a2 



f{v 



2vr ^l + Aa{a + v) 

The PF has a weak singularity at f = — l/(4a) — a 
2 exp(-^) 



(13) 



(14) 



^ ^Jl + Aaia + v) 

The PF is shown in Fig. 3 for a several values of the parameter a. The left hand side panel 
shows the PF for a relatively large amplitudes a = 0.25, 0.18, 0.1, 0.05. For small amplitudes 
a better illustration is the difference of the PF and Gaussian PF shown in the right hand 
side panel (a = 0.02, 0.01, 0.005). Note that for small a 

fiv)^[l + ahiv)]fGiv) (15) 



In the paper by Luo & Schramm ( 1993) the second term is missed which is not important for small a. 
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10 S.Shandarin 
where 

h{v) = viv"^ - 3) (16) 

and fdv) is the Gaussian PF. The function 0.02h{v)fG{v) is shown by the solid line in the 
right hand side panel of Fig. 3. It almost merge with the exact PF shown by the dotted line. 
The agreement is even better for smaller values of a. 
The two and three point functions are respectively 

6(1,2) = e,(l,2) + 2a%2(i^2), 
^3(1, 2, 3) = 2a [^^(1, 2)^,(1, 3) +^,(1, 2)^,(2, 3) +e,(l, 3)^,(2, 3)] 

+ 8a=^e.(l, 2)^,(1, 3)^,(2, 3). (17) 

In the case of quadratic mapping (|) the three-point function detects non-Gaussianity in 
the linear order in a. In this case all even order functions vanish in linear order of a due to 
symmetry of the mapping. 

The global Minkowski functionals of the field can be also readily obtained 



C{v) = Cg —j= \+Cg 



^l + 4:a{a + v)- 


-1 


2^/2a 




^l + 'ia{a + v) - 


- 1 



2v^, 



a 



where Cg{u) and Gg{u) are the Gaussian global Minkowski functionals (eq.0). The x^ field 
with one degree of freedom is obviously a particular limiting case of this model and has been 
discussed in great detail QTomita 1990| ; [Worsley 1994| ; [Coles fc Barrow 1987| ; pchmalzing 
19991 INovikov, Schmalzing fc MukhanovlOOOl) . 



Figure 4 (left hand side panels) shows global Minkowski functionals as a function of the 
level for four values of the amplitude a = 0.25, 0.18, 0.1, 0.05. Combining the equation for 
the CPF (|T2|) and recalling that A{v) = F{v) one can also plot C and G as functions of A 
(Fig. 4, right hand side panels). The Gaussian curves are shown by solid lines in all panels. 
One easily sees that the non-Gaussian signature is much stronger in the left hand side panels. 
This is due to leaking of some non- Gaussian information into the C{v) and G{v) curves from 
the PF. The curves C{A) and G{A) show only the non-Gaussianity that is absent in PF 
and therefore show nothing when such information is absent. This is why one can see only 
two non-Gaussian curves corresponding to a = 0.25 and 0.18 in the right hand side panels. 
Although all four non-Gaussian curves are plotted the curves corresponding to a = 0.1 and 
0.05 are merged with the Gaussian curves and are not seen. The curves corresponding to 
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Morphological non-Gaussianity 11 

a = 0.1 and 0.05 are clearly seen in the left hand side panels although, as shown below, 
the morphology of these fields is practically Gaussian or more exactly the non-Gaussian 
morphology is completely absent in maps with Npi^ < 10^. Another illustration of this effect 
is provided by the comparison of the C{v) and G{v) curves at t; > 1 with the C{A) and G{A) 
at y4 < 0.16. The mapping (eq. |]) is monotonic and thus does not change the morphology of 
the field at small amplitudes in this range of v or A. All the curves merge in the right hand 
side panels manifesting the similarity of the morphology to the Gaussian one, but they are 
clearly distinct in the left hand side panels due to leaking of some non-Gaussianity from the 
PF. 

As a whole the quadratic model (|^) is not a monotonic mapping: v is monotonically 
increasing at at m > —l/2a and monotonically decreasing at -u > —l/2a. This results in two 



terms in eq. |I2] and eq. |1^ as well as in the deviation of the morphology from the Gaussian 
one because the regions belonging to the excursion set v > Vc belong to different excursion 
sets in the parent field u: Si [u > 'Ui(fc)] and 5*2 [u < U2{vc)] (eq. [TT| ). 

Consider now the limiting case of small a. For small a the decreasing branch of the 
mapping is shifted to large negative u and thus greatly reducing the second terms in eq. (O 



and |18]) . Therefore 

A{v) ^ - 



^ l + 4:a{a + v) - 1 
1 — erf 



2^2. 



a 



C{v) ^ Cg 



2^2: 



a 



'1 +4a(a + v) -1\ 

and the mapping becomes almost monotonic. 

In cosmology one is usually interested in detecting the possibly smallest non-Gaussianity. 
For small a the quadratic transformation can be practically monotonic for a map of a finite 
size. Let us consider a map with Npix pixels and calculate the critical value ac corresponding 
to the value of a when the double valued character of transformation (P) is observed on 
average at one pixel only. The minimum of the parabola (^ is reached at Umm = — l/2a. 
Therefore the probability that a map of the size Npix contains on average one pixel where 
the monotonicity of (^ is broken is 

1 - Fciumin) = 1 - FGi-l/2a,) = -^. (20) 
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12 S.Shandarin 

The solution for Npi^ has a simple analytic form for a^ ^ 1 

iV«x ~ a/-— exp( — -). (21) 

Critical amplitude etc is plotted as a function of logiQ^Npi^) in Fig. 5 with the points cor- 
responding to the best resolutions of the COBE, MAP and Plank experiments. One easily 
sees that only relatively large a (roughly a > ac ^ 0.09) may result in the deviation of 
the quadratic model (|^) from monotonicity for even very large maps. It means that the full 
information about the non-Gaussianity of the quadratic model can be obtained from a one 
point function only if a < 0.09. Nothing can be gained by applying more complex n-point 
statistics. 

5 TRIVIAL NON-GAUSSIANITY 

From the examples of Sec. 4 one can make a further step to a little more general case. 
Consider a monotonic mapping of a Gaussian field u ^ v. v = 0('u) with dv/du > (the 
case with dv/du < is similar). This transformation obviously affect only the PF 

F{v) = FG{r\v)) (22) 

where Fq is the Gaussian cumulative probability function and (j)^'^{v) is the inverse of the 
function 0('u). The PF can be found by the differentiation of the above expression. 

The shapes of the level contours do not change in the resulting field because the contours 
of constant v coincide with the corresponding contours of constant u. The non-Gaussian 
fields obtained by this kind of mapping can be called trivial since all its non-Gaussianity is 
described by the one-point function (e.g. PF). 

The cubic as well as quadratic model with roughly a < 0.09 (in maps with Npi^ < 10^) 
are examples of trivial non-Gaussian fields. Another example of a trivial non-Gaussian field 
is given by the exponential mapping 

V = exp(u) (23) 

where m is a Gaussian field. The resulting field has the lognormal PF and exactly the 
Gaussian morphology. 

If f = <f){u) is not a monotonic function then u = (j)^^{v) is a multiple valued function. 
As a result some regions of the excursion set with u > Uc may not satisfy the condition 
v{u) > Vc{uc) and vice versa. This case is a little more complex but also can be fully treated 
analytically in terms of the Minkowski functionals . 
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Morphological non-Gaussianity 13 

One point function PF (or CPF) is in many respects the simplest statistics requiring 
only 0{Npix) operations. The other statistics (e.g. n-point correlation functions) share some 
information with the PF. The morphological statistics (Minkowski functionals ) are easy to 
isolate from the PF by parameterizing them by the area of the excursion set. We are not 
aware of the general method of isolation the morphological non-Gaussianity in other non- 
Gaussian statistics. However, there is a simple technique allowing to assess the sensitivity 
of any statistics to the morphological non-Gaussianity. 

Let us relabel the levels of the field f — > m according to the equation 

F{v) = Fg{u) (24) 

where F{v) and Fg{u) are the CPF of the original and Gaussian fields respectively. The new 
field has the Gaussian PF by design but restores the morphology of the original field because 
the contour lines have been only relabeled but not distorted. Then computing any non- 
Gaussian statistics for the new field one can be sure that it depends only on the information 
that is absent in the PF. For instance, if the original field was trivial non-Gaussian then 
after this transformation all non-Gaussian statistics must vanish. This transformation of the 
field can be called the gaussianization of the field. It was used for the purpose of recovery of 
the primordial fluctuations from the present day galaxy distribution ( [Weinberg 1992| ) and 



recovery of the power spectrum from QSO Lya forest spectra ( Croft et al. 1999| ). 



It is worth stressing that the gaussianization neither simply smoothes the original field 
no re-scales the spectrum as in the transformation suggested by Wu (| 2001 ). It is local 



but in general highly nonlinear and its only purpose is to isolate non-Gaussian information 
which is complimentary to the PF ( [Shandarin et al. 2001 ). It certainly affect all n-point 



functions and n-point spectra including the spectrum of the field. It also changes the phases 
and therefore may be useful in studies similar to that by Coles & Chiang (| 2000[ ). 



6 THE AMPLITUDE OF THE QUADRATIC MODEL 

Assuming that detecting the possibly smallest non-Gaussianity is the goal we consider a 
few estimators of the amplitude of the quadratic non-Gaussianity a ignoring minute effects 
on morphology. It means that a is to be less than Uc of eq. ^ (see also Fig. 5). Luo and 
Schramm (1993) derived the skewness in the quadratic model. To linear approximation in 
a it is 

© 0000 RAS, MNRAS 000, 000-000 



14 S.Shandarin 

S ^ 5i!iZ^ =< ,3 >^ g^_ (25) 

(J 

The last equality is because w = and a^ = 1. Thus, measuring the skewness one can 
estimate the amplitude of the quadratic component. One can also easily obtain the estimates 
of a from higher central moments 

C5 ^ 60a, C7 ^ 630a (26) 

where q are the corresponding central moments q =< f * > since < f >= and < f^ >= 1. 
In principle, any of odd moments can be used for the estimate of the amplitude a. 

Since the higher moments depend stronger on the tails of the PF one might think that 
the most accurate is the estimate based on the skewness. As we shall see this trend is very 
week at least in the case of the lowest three moments (n=3, 5 and 7). 

Using the asymptotic form of f{v) (eq. |14D one also can derive a using the least square 
estimator. Suppose there are N measurements of v binned in N^ equal bins, Sv, {rii, i = 
1, ■ ■ ■ , Nh). Then minimizing the x^ 



N, 



x^ = E 



i=l 



Hi -n{a,Vi)) 



(Ji 



(27) 



one can obtain a. Simple calculations result in 

Considering the propagation of errors one can also estimate the standard deviation of a 



Assuming rij ^ 1 one can simplify eq.29 



(N^ \-l/2 / oo \-l/2 

<y. ^ iV-V2 \y^fG{v,)h\v.:)bv\ ^ N-'l^ j fG{v)h\v)dv\ = {QN)-'I\ (30) 

This can be used for the estimate of the statistical limit on the value of a on a grid of size 
Npix- Requiring that |a| > 2aa (2(j detection) one obtains 

\a\ > {m,,j2Y^'\ (31) 

Thus, for the map of the size of Npi^ ~ 10^ the value a must be greater than amin ~ 10~^ in 



order to be detected by this method. For an actual CMB experiment eq. |30| is good but still 
an approximation. The effects of the radiation transfer function are needed to be included 
in a realistic application to CMB map. 
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7 MONTE CARLO SIMULATIONS 

In order to check the above analytic calculations and study the other effects we use simple 
Monte Carlo simulations. First, a Gaussian field u {<u >= 0, < u"^ >= 1) is generated in a 
square of the size Npi^ = Ni x Ni. We used the fiat spectrum smoothed at ks = Ni in units 
of the fundamental wavelength, P oc exp(— A;^/A;^). Then, the field v = u + a(-u^ — 1) was 
obtained and the amplitude a was estimated by four methods described above: ai from the 
least square estimator (eq. ^), a^ from the skewness (eq. ^^, and a^ and ay from the fifth 
and seventh moments respectively (eq. |2^ ). 

In addition, the mean of four was also computed: am = (ai + as + as + a7)/4. Figures 
6 and 7 show the probability functions P{<y) measured in 400 realizations of the 256^ and 
1024^ maps. Four panels show results for four different amplitudes illustrating the reliability 
of the detection of the non-Gaussianity and measurement of the value of a. Five curves 
in each panel show the distribution functions P(ai), P^a^), P^a^), P{a7), and P{am)- It 
is a little surprising that the curves heavily overlap which demonstrates that there is no 
much difference in the distributions of a obtained by different methods although the mean 
of four is always has a little shorter tails (solid lines) while the estimate based on the highest 
moment Cy often spreads a little wider than others (dotted lines). Although all a, are not 
independent they are the most sensitive to different parts of the PF. The a-i is obviously 
sensitive to the higher values of the PF than others and thus is probably the least reliable. 

Comparing the theoretical estimates of the statistical limits on the measurement of the 
quadratic non-Gaussianity (0.003 for the 256^ and 0.0008 for 1024^ map eq. ^0|) with Fig. 6 
and 7 we conclude that eq. ^ predicts the errors in measuring a quite accurately. 

In these simulations the field v obviously has a scale that corresponds roughly to the pixel 
size. If the scale of the field was below the resolution of the map the field would be smoothed. 
Smoothing obviously erases non-Gaussianity because of the central limit theorem. We model 
the resolution effect by smoothing the original non-Gaussian field with the Gaussian filter 

W{k,hm)=exp{-^). (32) 

sm 

Real CMB maps are likely to be smoothed which obviously reduces the ability of measur- 
ing a. Figure 8 illustrates this effect. The quadratic non-Gaussianity can be easily detected 
and the amplitude a = 0.01 can be roughly estimated from a 256^ map if the scale of 
non-Gaussianity equals the pixel size. In the top left panel the distribution function of a 
measured in a thousand maps peaks at about right value a = 0.01. The width of the distri- 
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bution function definitely allows the detection but the accuracy of the measurement is not 
very high. If it was observed with four times lower resolution corresponding to the smooth- 
ing scale ksm = 64 in the units of the fundamental wavenumber the measured amplitude 
would be on average a half of the original one. At the resolution 8 times of the original 
scale of the non-Gaussianity the detection becomes impossible. Of course, this depends on 
the amplitude. If the amplitude was 2.5 times higher {a = 0.025) then the detection of the 
non-Gaussianity would still possible even at the resolution 8 times of the original scale as 
illustrated by by the right top panel of Fig. 9 however the measurement of the amplitude 
would not be possible even it was ten times higher a = 0.1. Generally, smoothing reduces 
non-Gaussianity and thus makes its measurement more difficult. 

8 SUMMARY 

The simplest tests for the primordial non-Gaussianity are based on studies of the one-point 
functions: the PF or CPF and in particular the moments of the PF (skewness, kurtosis, 
etc. ) . However, the Gaussian form of the PF does not guarantee the Gaussianity of the field. 
The hierarchy of the n-point functions or n-spectra or other tests are supposed to analyze 
additional information that is not contained in the PF. It is natural to call it morphological 
information because it is determined by the pattern of the contour lines only. 

However, the sensitivity of n-point functions as well as many other statistics to the 
morphological information is unknown. In fact the n-point functions and other statistics 
use a great deal of the information stored in the PF. A simple example is provided by the 
analysis of the quadratic model (eq.|]). 

At small amplitudes (practically at a < 0.1 Q]) the field becomes trivial non-Gaussian 
which means that the full information about the non-Gaussianity is already contained in 
the PF (eq. |15D. 

This means that nonzero three-point function if applied to this model (eq. p!7|) uses 
exactly same statistical information as the PF. The morphological statistics in the form of 
the Minkowski functionals of the excursion set {A, C and G) and percolating region {Ap, Cp 
and Gp) obviously recycle the same information as the PF. However, in the last case this 
information can be easily isolated in A which is the CPF of the field. Other morphological 

'I One can see that even the value ac = 0.08 requires a map with Npi^ > 10® (Fig- 5) and eq. Gil shows that ac is extremely 
weak function of Npi^: . 
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statistics {C, G Ap, Gp and Gp) do not show any non-Gaussian signal if they are expressed as 
functions of A (Fig. 4). In general, these morphological statistics use the different statistical 
information than the one-point functions. 

Expressing the morphological parameters as functions of A instead of the level u is 
equivalent to the gaussianization of the field. For other (non morphological) statistics the 
gaussianization of the field v can be achieved by relabeling the levels of the field v ^ u 



according to eq. |2^. The new field u has the Gaussian PF by design but the geometry and 
topology of the original field. If any statistic sensitive to non-Gaussianity is applied to the 
new field u will characterize the geometry and topology of the original field v but free from 
the effects of the non-Gaussian PF of the original field. Performing this test allows to assess 
the level of independence of any statistic from the one-point functions. 

In line of the above reasoning we showed that the quadratic model results in a trivial 
non-Gaussian field at small a. Using this fact we estimated the statistical limit of detecting 
the quadratic model in the ideal situation when the scale of the field is resolved. This limit is 
roughly asi ~ Np^J if the parent field u is normalized to unity as in eq. ^ or a'^icr' ^ Np^J 
if the rms of the parent field is a'. The Monte Carlo simulations on 256^ and 1024^ grids 
confirm this estimate. 

We measured the amplitude of the quadratic non-Gaussianity a using four different 
estimators of a: ai from the least square fitting of the PF (eq. ^8]) and from three lowest 
order odd central moments 5 = C3, C5 and cj 

S C5 Cj , . 

"^=6' "^ = 60' "^=630- ('') 

In addition, we also computed the mean of the four, am = (ai + as + 05 + a7)/4. Quite 
surprisingly, all four estimators ai, 0:3, 0:5, 07 performed very similarly in all tests. This fea- 
ture in combination with the Gaussian character of morphological statistics can be used for 
distinguishing the quadratic non-Gaussianity from other types of non-Gaussianity. 

We have illustrated how the resolution and smoothing can affect the detection and mea- 
surement of the quadratic non-Gaussianity. In order to be detected the scale of the original 
non-Gaussian field must be relatively close to the resolution of the map. However, quantify- 
ing this effect requires the study of three-dimensional parameter space (the size of the map, 
Npix, the amplitude of the quadratic non-Gaussianity, a, and the resolution or smoothing 
scale, ksm) and is beyond the scope of this paper. 
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Figure 1. Global Minkowski functionals of a Gaussian field. Three panels on the left hand-side show three global Minkowski 
functionals: the total fraction of the area in the excursion set, A, contour length per unit area, C, and genus per unit area, G, 
as a function of the level, u. The panels on the right hand-side show (from top to bottom) the level, n, contour length, C, and 
genus G as a function of the total area A. 





Figure 2. Percolation curves of Gaussian fields. Three panels on the left hand-side show: the area, Ap, contour length, Cp and 
genus, Gp of the largest region as a function of A in the same units as in Fig.l. The right hand side panels show the ratios of 
the Minkowski functionals of the percolating region to the corresponding global Minkowski functionals . The error bars show 
Iff errors of the Monte Carlo simulations 
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Figure 3. The probability functions in the quadratic model for five amplitudes: a = 0.25,0.18,0.1 and 0.05 corresponding to 
the dot-dashed, long dashed, short dashed, and dotted lines are shown in the left hand side panel. The solid line is the Gaussian 
PF. The difference between the PF of the quadratic model and the Gaussian PF for three amplitudes : a = 0.02, 0.01 and 
0.005 corresponding to dotted, long dashed and short dashed lines is shown if the right had side panel. The solid line shows 
the analytic linear approximation (eq. hStj with the amplitude a = 0.02. 



0.6 




Figure 4. The Global Minkowski functionals C and G for the quadratic model parameterized by the level v are shown in 
the left hand side panel, the parameterization by the area of the excursion set A is shown in the right hand side panels. The 
amplitudes are as in the left hand side panel of Fig. 3: a = 0.25,0.18,0.1 and 0.05 corresponding to the dot-dashed, long 
dashed, short dashed, and dotted lines. The solid lines show the Gaussian Minkowski functionals . 
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Figure 5. The triviality limit for the quadratic model ac = cic{Npix). If the amplitude a < Oc then the quadratic model is 
fully described by the PF only on a grid of the size Npi^ . 




300 I I I I 



0.005 0.01 



0.005 0.01 



Figure 6. The probability functions of the amplitudes measured in the Monte Carlo simulations of the 256^ maps. The 
amplitudes oi (obtaind from eq. KSl) , «3 (eq. ba), 05,07 (eq. Eq) and Om = (ai +03 + 05 + 07)/4 are shown by the dot-dashed, 
long dashed, short dashed, dotted, and solid lines respectively. The vertical dotted lines mark o = and the model value am 
in the simulations. 
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Figure 7. The probability functions of the amplitudes measured in the Monte Carlo simulations of the 1024 map. Notations 
are as in Fig. a 
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Figure 8. The effect of resolution in the 256 maps. Four panels show the PFs of the measured amplitudes when the initial 
quadratic field was unsmoothed (left top), filtered with the Gaussian window with the Nyquist frequency ksm = 128 (left 
bottom), with fc^m = 64 (right top), and with ksm = 32 (right bottom). Notations are as in Fig. a 
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Figure 9. The effect of the amplitude in the 256 maps. Four panels show the PFs of the measured amplitudes if the initial 
quadratic field smoothed with fcsm = 32 had the amplitudes a = 0.1 (left top), a = 0.05 (left bottom), a = 0.025 (right top), 
and a = 0.01 (right bottom). Notations are as in Fig. H 
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